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A simple, yet subtle 



of two-body decay kinematics 
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We study the two-body decay of a mother particle into a massless daughter. We further assume 
that the mother particle is wnpolarized and has a generic boost distribution in the laboratory frame. 
In this case, we show analytically that the laboratory frame energy distribution of the massless decay 
product has a peak, whose location is identical to the (fixed) energy of that particle in the rest frame 
of the corresponding mother particle. Given its simplicity and "invariance" under variations of 
the boost distribution of the mother particle, our finding should be useful for the determination 
of masses of mother particles. In particular, we anticipate that such a procedure will then not 
require a full reconstruction of this two-body decay chain (or for that matter, information about 
the rest of the event). With this eventual goal in mind, we make a proposal for extracting the 
peak position by fitting the data to a well-motivated analytic function describing the shape of such 
energy distribution. This fitting function is then tested on the theoretical prediction for top quark 
pair production and its decay and it is found to be quite successful in this regard. As a proof of 
principle of the usefulness of our observation, we apply it for measuring the mass of the top quark 
at the LHC, using simulated data and including experimental effects. 

PACS numbers: ll.80.Cr 

It is very well-known that in the rest frame of a mother 
particle undergoing a two-body decay, the energy of each 
of the daughter particles is fixed in terms of mother and 
the daughter particle masses (see, for example, Ref. [1] 
and the references therein). Turning this fact around, 
we can determine the mass of the mother particle if we 
can measure these rest-frame energies of the daughter 
particles. 

However, often the mother particle is produced in the 
laboratory with a boost, that too with a magnitude and 
direction which is (a priori) not known. Moreover, the 
boost of mother particles produced at hadron colliders 
is different in each event. Such a boost distribution de- 
pends on the production mechanism of the particle and 
of the structure functions of the hadrons in the initial 
state of the collision and is thus a complicated function. 
In turn, the fact that the mother has a different boost 
in each event implies that when we consider the observed 
energy of the two-body decay product in the laboratory 
frame, we get a distribution in it. Thus it seems like the 
information that was encoded in the rest frame energy 
is lost and we are prevented from extracting (at least at 
an easily tractable level) the mass of the mother particle 
along the lines described above. 

We show that, remarkably, if one of the daughter parti- 
cles from the two-body decay is massless and the mother 
is unpolarized, then such is not the case. Specifically, 
in this case, we demonstrate that the distribution of the 
daughter particle's energy in the laboratory frame has a 
peak precisely at its corresponding rest-frame energy. 

This result is interesting per se. Furthermore, we ex- 
pect that it will lead to formulation of new methods for 
mass measurements. Obviously, for this purpose, we need 
to be able to determine the location of this peak accu- 
rately from the observed energy distribution of the mass- 



less daughter. To this end, we propose and motivate an 
analytic function that can be used to fit the data on the 
energy distribution and thus extract the peak position. 
We show that this function is a suitable one using the 
top quark decay, t — > W~b, as a test case, namely, it fits 
very well the theory prediction for energy spectrum of 
the resulting 6-jets. Simulating a realistic experimental 
situation, we then show that we can extract the value of 
the top mass from the position of the peak in the 6-jet 
energy distribution along with the well-measured mass of 
the W boson. 



Let us consider the decay of a heavy particle B of mass 
nig, i.e., B — > Aa where a is a massless visible particle. 
For the subsequent arguments, the properties of the par- 
ticle A (other than its mass which we denote by tua) are 
irrelevant. In the rest frame of particle B, the energy of 
the particle a is simply given by: 



E* = 



2m i 



(1) 



Here and henceforth the starred quantity denotes that it 
is measured in the rest frame of particle B, i.e., mother 
particle. If the mother particle (originally at rest) is 
boosted by a Lorentz factor 7 = 1/yl — (3 2 in going to 
the laboratory frame, then the energy of particle a seen 
in the laboratory frame is 



E = £*7(1 + /3cos0*) 



(2) 



where 9* defines the direction of emission of particle a in 
the rest frame of B with respect to the boost direction 
/? of the mother B in the laboratory frame. Due to our 
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assumption of the mother being not polarized 1 the distri- 
bution of cos 9* is flat. This implies that the distribution 
of E is flat as well. More precisely, since cos 9* runs from 
— 1 to +1, for any fixed 7 the shape of the distribution 
of E is a simple "rectangle" spanning the range 2 



E 

E* 



(7 - vV - 1) . (7 + W - 1) 



(3) 



where x defines the dimensionless energy variable of the 
visible particle in the laboratory frame normalized by 
its rest-frame energy (here, we have used the relation- 
ship, 7/3 = \J^j 2 — 1). A few crucial observations arc 
in order. First, the lower (upper) bound of eq. (3) is 
smaller (larger) than 1, which implies that each rectan- 
gle for a fixed, but arbitrary 7 always contains E* . In 
turn, this implies that, for any distribution of the boost 
of the mother particle J5, the energy distribution of the 
daughter a always covers the energy E* . Remarkably E* 
is the only value of the energy to enjoy such a property, 
as long as the distribution of mother particle boost is 
non- vanishing in a small region around 7=1. Thus, the 
peak of the energy distribution of the particle a is unam- 
biguously located at E = E* . In fact, this argument goes 
through even for a massive daughter, provided we restrict 
boosts of the mother particle to 7 < (27* 2 — l) , where 
7* denotes the boost of the daughter in the rest frame 
of the mother. Secondly, such rectangles are asymmetric 
with respect to the point x = 1 (where E — E*), i.e., 
the upper bound is farther from x — 1 than the lower 
bound. Thus, the energy distribution of the particle a 
has a longer tail toward high energy w.r.t such a peak. 

More formally, one can derive the (normalized) differ- 
ential decay width in x from the fact that cos9* is flat 
as mentioned above. The result for a fixed 7 is 



1 dr 

r dx 



e (x - 7 + - 1) © (-x + 7 + vV - 1) 



fixed 7 



2VT 2 



(4) 



where Q(x) is the usual Heaviside step function, and the 
two step functions here merely define the allowed range 
of x. Next, consider a probablility distribution of boosts 
of the mother, given by 5(7). A given energy of daughter 
in laboratory frame (x) can actually result from a specific 
range of values of the mother boost (7), as per eq. (4). 
So, we have to superpose these contributions, weighted 
by the boost distribution, giving: 



/(*) 



1 dT 

r dx 



G?7 



5(7) 



2^f 



(5) 



The lower end in the integral here was derived from the 
solution to the following equation for 7: 



x = 7 ± \/j 2 — 1 . 



(6) 



Here the positive (negative) signature is relevant for 
x > 1 (x < 1). From eq. (5) we can also compute first 
derivative of f(x), that is 



sgn(l ~ x) 
2x 



1 



^2 [ X+ x 



(7) 



We assume that 5(7) does not vanish for any finite 
non-zero value of 7. This is rather typical for particles 
produced at colliders. In what follows we show that this 
is sufficient to guarantee that there is a peak at E* . We 
consider the two possibilities for g(l). Namely, if it van- 
ishes, then f'(x = 1) oc g(l) = 0, and the distribution 
has its unique extremum at E — E* , following from our 
assumption on .9(7). If g(l) ^ 0, then f'(x) flips its sign 
at x = 1 so that the energy distribution has a cusp at 
E = E* . The function f{x) is positive and vanishes for 
both x — > and x — > 00, therefore the point E — E* is 
necessarily the peak of the distribution in both cases. 

The presence of the peak is therefore demonstrated 
(both heuristically and formally). In particular we have 
shown that the peak position does not depend on the 
actual boost distribution of the mother particle, as long 
as this distribution is non-trivial, i.e., 5(7) is not a <5- 
function. 

As advertised at the beginning, our finding can be uti- 
lized to measure a combination of and mg given by 
E*, which means that tub can be measured if itia is 
known, or vice versa. For this purpose, we are required 
to extract the location of the peak accurately from data. 
Clearly, having a theoretical prediction from first princi- 
ples for the shape of f(x) is very hard because the boost 
distribution 3(7) depends on the structure functions and 
on the actual decay vertex of the mother particle, and 
hence in principle it depends on the process, the mass of 
the mother particle, and the collider. Nevertheless, we 
can still exploit some functional properties of the generic 
f(x) which are listed below: 

• / is a function with an argument of \ (x + ^J, i.e., 
it is even under x •<-> — , 

• / is maximized at x = 1, 

• / vanishes as x approaches or 00, 

• / becomes a (^-function in some limiting configura- 
tion. 



1 For example, the mother can be a scalar particle or even if it 
has spin, then its production process might not distinguish the 
various polarizations. 

2 This result is quite well-known and was used for a measurement 
of the W mass at lepton colliders [2]. 



The first property follows from eq. (5), and the second 
from eq. (7) and the argument thereafter. The third one 
is also manifest from eq. (5) since those two limits lead to 
a trivial definite integral from 00 to 00. Finally, the last 
one reflects the fact that no boost returns a fixed value of 
energy given in eq. (1), i.e., (S-functionlike distribution. 
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FIG. 1: The orange dots are the theory prediction for da/dEb 
in the process eq. (10) computed with MadGraph5 at LHC with 
yfs=7 TeV. The purple line is the best fit of our ansatzeq. (8). 
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FIG. 2: An instance of the result of the fit on the energy 
distribution of the 6-jets in a pseudo-experiment. For the fit 
we binned data in bins of 4 GeV. Only the blue data points 
are used in the fit, which correspond to using only the part 
of the spectrum from 30 to 150 GeV. 



Being aware of the constraints given above, we propose 
the following "simple" function as an ansatz for f(x): 



f(x) = K x ^exp 



P f , 1 

-- \ x H — 
2 \ x 



(8) 



where p is a parameter which encodes the width of the 
peak and the normalization factor K\(p) is a modified 
Bessel function of the second kind of order 1. One can 
easily prove that the proposed ansatz can be reduced to 
a (5-function for any sufficiently large p using the asymp- 
totic behavior of Ki(p) such that 



p— >oo 



lp 



a 1 - 

,P 



(9) 



Such a behaviour is required in order to reproduce the 
limit of no boost of the mother. Finally, we can show 
that the above ansatz does not have a cusp (at E*) so 
that it is more suitable for the case of g(l) — such as 
pair-production of mothers 3 . 

In order to test the goodness of the ansatz given in 
eq. (8) we use it to fit a theoretical prediction for the 
distribution of 6-jet energy in top decay. The bottom 
quark is not massless; it is nonetheless highly boosted in 
the rest frame of top quark, namely, 7* ~ 15. Based on 
our earlier discussion of the massive case, our argument 
for the peak in 6-jet energy being at E* is invalidated 
for boosts of the top quark which are so large (7 > 500) 
as to have a negligible probability. Hence, we expect the 
peak to be very close to E* . Similarly we expect that the 
first of the functional properties of the energy spectrum 
eq. (5) will be only negligibly violated by the non-zero 



3 We defer the study of the other case of g(l) 7^ 0, for example, 
the single production of the mother, for future work. 



mass of the bottom quark. This justifies the use of the 
ansatz eq. (8) to fit 6-jet energy spectrum. 

Specifically, we study a sample of fully leptonic top 
decays from the process 



pp — > tt — > bbfi 



(10) 



at the Large Hardron Collider (LHC) with 7 TeV center- 
of-mass energy. To compute the theory prediction for 
the process eq. (10) we employ MadGraph5 1.4.2 [3] as a 
matrix element generator. We also take m top of 173 GeV 
and the patron distribution functions (PDFs) CTEQ6L1 [ ] 
evaluated with a renormalization and factorization scale 
varied depending on the kinematics of each event accord- 
ing to the default of MadGraph5. 

The result of the associated fit is exhibited in FIG. 1 
which shows a very good agreement between the the- 
ory prediction from MadGraph5 and the fitting function. 
To quantify the goodness of the ansatz with an objec- 
tive measure we compute both the Kolmogorov-Smirnov 
(KS) [5] and the x 2 value. The latter is computed tak- 
ing bin counts for a luminosity of 5/& -1 at LHC with 
y/s = 7 TeV assuming that the error on each bin count is 
gaussian. The result is \ 2 — 39.3 for 198 degrees of free- 
dom while the KS test statistic is 0.012. Neither the % 2 
nor the KS test has any particular probabilistic meaning 
attached to it, in fact they just serve the purpose of quan- 
tifying in an objective manner the agreement between the 
theory prediction and the ansatz. The \ 2 should be sen- 
sitive to deviations in regions of the distribution where 
there are more events, i.e., at the peak, while the KS test 
should be sensitive to the overall shape. Both tests indi- 
cate that our ansatz gives a very good fit to the theory 
curve. Although not shown in the figure, the agreement 
is very good also in the high-energy tail. We have inves- 
tigated the sensitivity of this result to the choice of the 
PDFs by repeating the same fit on the theory prediction 
obtained using the MRST2002NL0 PDFs set of Ref. [ ]. We 
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observe negligible differences w.r.t. the result obtained 
with CTEQ6L1. 

So far, we have found that the ansatz in eq. (8) is very 
good at reproducing the theory prediction. In fact, this 
success suggests that the ansatz may be used to measure 
the combination of masses in eq. (1) from experimental 
data. In order to investigate this possibility, we go back 
to the example of the top quark, namely, we would like 
to use the fitting function in order to extract the peak of 
the observed energy distribution of the 6-jet and measure 
the top quark mass by plugging this value and the well- 
known mass of the W boson into eq. (1). Before getting 
into details, we would like to mention that, we do not nec- 
essarily aim at getting a result for the value of mtop that 
is competitive with the current measurements. Rather 
we aim at finding what is the sensitivity of our method 
for measuring top quark mass in a realistic setup. In fact 
for a fair comparison it should remarked that the current 
measurements of mt op rely on rather complicated tools 
and often advocate templates for the distributions that 
require a detailed knowledge of the underlying dynamics 
of the top quark decay, i.e., the matrix element of the pro- 
cess. On the contrary, our method is extremely simple: 
it is based on pure kinematics and does not rely at all on 
detailed knowledge of the above-mentioned dynamics (as 
long as the top quark is produced unpolarized) . As such 
we can regard our study of the mass measurement of the 
top as a proof of principle that our method can be used to 
measure the mass of heavy particles, in particular, new 
physics particles. 

In lieu of actual experimental data we use a sam- 
ple of Monte Carlo (MC) simulated collision events. 
Namely, we further process the event sample generated 
by MadGraph5 1.4.2 to include the effects of showering 
and hadronization as described in PYTHIA 6.4 [7] with 
detector response simulated by Delphes 1.9 [8] and jets 
made with FastJet [9, 10] using the anti-kr algorithm 
[11] with the parameter choice R — 0.4. Furthermore, we 
impose cuts on the final state of the process eq. (10) fol- 
lowing the selections of Ref. [12] for the e/j, final state. We 
consider an ensemble of 100 pseudo-experiments, each of 
which is equivalent to of data from the LHC at 

^fs = 7 TeV. For each pseudo-experiment we perform a 
fit with our template eq. (8). From the extracted value 
of the peak of the distribution we get a measurement of 
m top . The distribution over the 100 pseudo-experiments 
of mtop and its la error are symmetric around the cen- 
tral values, and do not show special features. For a bin 
size of 4 GeV the average best-fit m top and the la error 
resulting from the fit are 

(m t0 p> = 173.1 ± 2.5 GeV (11) 

with a median of the reduced x 2 of our fit equal to 1.1. 
We have checked that the result of the fit is stable un- 
der changes of the range of energies taken in the fit and 
changes of the bin size. For illustration purpose we show 
the outcome of one of the pseudo-experiments in FIG. 2. 



The obtained value of m t0 p is rather good: the er- 
ror is small and the central value is compatible with the 
value from the standard methods (within the la error of 
the latter). Furthermore, the obtained x 2 is good. All 
this tells that the usefulness of our function eq. (8) is 
not spoiled by the selection criteria for top quark decay 
events nor by detector effects. Even though our assess- 
ment of the power of the technique does not take the 
background and the entire realm of detector effects into 
consideration, we regard it as rather encouraging for the 
determination of the masses of heavy particles at colliders 
based only on the kinematics of the decay. 

In conclusion, we have shown that for the two-body de- 
cay of an unpolarized boosted mother particle, the energy 
spectrum of a massless daughter in the laboratory frame 
encodes in a rather simple manner information about the 
masses involved in the decay. Specifically, we showed 
that, for a generic distribution of the boost of the mother 
particle, the location of the peak in the energy distribu- 
tion of the visible daughter is identical to its rest-frame 
energy. 

We also found a simple function that seems to repro- 
duce rather well the energy spectrum of the massless 
daughter. In particular, we assessed the goodness of this 
function by comparing it against the theory prediction 
for the energy distribution of fe-jets in the leptonic decay 
of top quark pairs. 

Furthermore we have studied how well our simple func- 
tion can be used to reproduce the distribution from MC 
simulated events, where the effect of soft QCD radiation 
and detector effects have been included. We found that 
our template function can be successfully used even in 
the presence of such effects. In fact, our finding can be 
used to extract the top mass in the LHC data with a 
rather good accuracy, well below 10%. Note that our 
method exhibits a remarkable twist in the mass mea- 
surement paradigm. Instead of using longitudinally or 
fully Lorentz invariant quantities for this purpose, we ex- 
tracted the mass of the top quark from a Lorentz- variant 
observable, i.e., energy of 6-jet. The crucial point is that 
even though the distribution of this quantity is depen- 
dent on the possible boosts of the mother particle, the 
location of the peak in it is "invariant" . 

While currently our technique might not compete with 
the standard (well-developed over more than a decade) 
methods to measure TOtop we remark that this observa- 
tion can be used to get a rather precise mass measure- 
ment of any heavy particle, in particular, new ones. The 
merit of our method is that it does not rely on any mea- 
surement of the other particle of the two-body decay so 
that we can extract some information about masses even 
if the latter is invisible, a case that would be not tractable 
with, for example, measurement of invariant mass only. 
It is also clear that our method and the traditional tech- 
niques for mass measurement are sensitive to different 
kind of detector effects. In general, we envisage that 
there will be a large degree of complementarity of our 
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method with more traditional ones. 

Finally, we emphasize that the proposed technique, de- 
spite being based on a fitting function, relies only on the 
minimal assumptions of absence of polarization and and 
the presence of a non-trivial boost distribution of the 
mother particle, i.e., it does not require any other prior 
knowledge about the underlying physics model governing 
the decay of the particle whose mass we want to measure. 
This suggests that our method will be especially suitable 
for the mass measurement of new particles to be discov- 
ered at the LHC, where we (a priori) would not know 
such details. 

In forthcoming work we shall describe a number of ap- 
plications of our finding here about the distributions of 
the energy of a daughter particle from a two-body decay. 
In particular, we shall present results on the measure- 
ment of the masses of specific new physics particles [13]. 
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